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Abstract 

o 

I In this paper we study the approximation of the distribution of Xt Hilbert-valued 

^SJ . stochastic process solution of a Unear parabolic stochastic partial differential equation 

written in an abstract form as 

O . 

O ■ dXt + AXt dt = Q^/^dWt, Xo = xeH, t e [0, T], 

. driven by a Gaussian space time noise whose covariance operator Q is given. We assume 

that A~°' is a finite trace operator for some a > and that Q is bounded from H into 
. D{A^) for some (3 > 0. It is not required to be nuclear or to commute with A. 

' The discretization is achieved thanks to finite element methods in space (parameter 

. h > 0) and implicit Euler schemes in time (parameter At = T/N). We define a discrete 

^ ' solution and for suitable functions ip defined on _ff , we show that 

|Ev3(Xf ) -E(^(Xt)| - ©(/i^"^ + Ai-^) 

where 7 < 1 — a + /3. Let us note that as in the finite dimensional case the rate of 
convergence is twice the one for pathwise approximations. 
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1 Introduction 



In this article, we study the convergence of the distributions of numerical approximations of 
the solutions of the solutions of a large class of linear parabolic stochastic partial differential 
\2{ ' equations. The numerical analysis of stochastic partial differential equations has been 

5^ . recently the subject of many articles. (See among others [I], [8], [H], [l2], [13], [H], [E], 

[l6],[II],[l8],[ig,[M!,[l5],[28],[29],[35],[M^ [Mj)- In all these papers, the aim is to give 
estimate on the strong order of convergence for a numerical scheme. In other words, on the 
order of pathwise convergence. It is well known that in the case of stochastic differential 
equations in finite dimension, the so-called weak order is much better than the strong order. 
The weak order is the order of convergence of the law of approximations to the true solution. 
For instance, the standard Euler scheme is of strong order 1 /2 for the approximation of a 
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stochastic differential equation while the weak order is 1. A basic tool to study the weak 
order is the Kolmogorov equation associated to the stochastic equation (see [22], [26], [27] 

In infinite dimension, this problem has been studied in much less articles. In [3j, the case of 
a stochastic delay equations is studied. To our knowledge, only [9], [20] consider stochastic 
partial differential equations. In [9] the nonlinear Schrodinger equations is considered. In 
the present article, we consider the case of the full discretization of a parabolic equation. 
We restrict our attention to a linear equation with additive noise which contains several 
difficulties. The general case of semilinear equations with state dependent noise present 
further difficulties and will be treated in a forthcoming article. This case is treated in 
|20j but there only finite dimensional functional of the solution are used and the finite 
dimensional method can be used. 

Note that there are essential differences between the equations treated in [3j and [9j. Indeed, 
no spatial difference operator appear in a delay equation. In the case of the Schrodinger 
equation the linear evolution defines a group and it is possible to get rid of the differential 
operator by inverting the group. Furthermore, in |9j, the data are assumed to be very 
smooth. 

In this article, we get rid of the differential operator by a similar trick as in [9]. However, 
since the linear evolution operator is not invertible, this introduces extra difficulties. More- 
over, we consider a full discretization using an implicit Euler scheme and finite elements for 
the spatial discretization. We give estimate of the weak order of convergence with minimal 
regularity assumptions on the data. In fact we show that as in the finite dimensional case 
the weak order is twice the strong order of convergence, both in time and space. 

In dimension d = 1,2, 3, let us consider the following stochastic partial differential equation: 
(1.1) ^!f^_Au(x,t) = 77(x,t), 

where x G O, a bounded open set of W^, and t g]0, T], with Dirichlet boundary conditions 

dr] 

and initial data and n = — with n denotes a real valued Gaussian process. It is convenient 

ot 

to use an abstract framework to describe the noise more precisely. Let be a cylindrical 
Wiener process on L^(C'), in other words is the space time white noise. Equivalently, 
given an orthonormal basis of L^(C'), W has the following expansion 

where (/3j)je?N is a family of independent standard brownian motions (See section 2.4 be- 
low). We consider noises of the form r/(t) = Q^/'^W{t) where Q is a non negative symmetric 
bounded linear operator on L'^{0). For exampl^, given a function q defined on O, we can 
take 

rj{x,t) = / q{x - y)W{y,t)dy, 
Jo 

Then the process ij has the following correlation function: 

Kr]{x,t)ri{y, s) = c{x — y){t A s) with c(r) = / q{z + r)q{z) dz. 

Jo 



^The equations describing this example are formal, it is not difficult to give a rigourous meaning. This 
is not important in our context. Note that q does not need to be a function and a distribution is allowed. 
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The operator Q is then given by Qf{x) = c{x — y)f{y) dy. Note that if the q is the Dirac 
mass at 0, r] = W and Q = I. 

Let us also set A = -A, D{A) = H^{0) n Hl{0) and H = L^{0). Then A : D{A) H 
can be seen as an unbounded operator on H with domain D{A). Our main assumption 
concerning Q is that A'^^'^Q defines a bounded operator on L'^(0) with a > —1/2 if d = 1, 
0" > if d = 2 and cr > 1/2 if d = 3. In the example above this amounts to require that 
(— A)'^/^(7 E L^{0). It is well known that theses conditions are sufficient to ensure the 
existence of continuous solutions of p.ip . 

If we write u{t) = u{-,t) seen as a -valued stochastic process then (ll.ip can be rewritten 
under the abstract Ito form 

(1.2) du{t) + Au{t)dt = Q^/^dW{t). 

In this article, we consider such an abstract equation and study the approximation of the 
law of the solutions of ()1.2p by means of finite elements of the distribution of n in ff. Let 
{^h}h>o be a family of finite dimensional subspaces of D{A^^'^). Let > 1 an integer and 
At = T/N. The numerical scheme is given by 

(1.3) {ul+'-ulvh) + At{Aul+',Vh) = VA't{Q'/\^+\vj,), 

for any Vh G V/j, where ^/Atx""^^ = W{{n + l)At) — W{nAt) is the noise increment 
and where (•, •) is the inner product of H. The unknown is approximated at time nAt, 
<n < N by ul eVh. In ([O]), we have used the notation ul'^'^ for Oul^^ + (1 - 0)ul for 
some 6 G [0, 1]. We prove that an error estimate of the following form 

|E((/?(n(nAt)) - E{ip{u';^))\ < c{h^^< + Af) 

for any function Lp which is and bounded on H. With the above notation, 7 is required 
to be strictly less than 1 — d/2 + a/2. This is exactly twice the strong order (see [33], [35]). 
If d = 1 and a = 0, the condition is 7 < 1/2 and we obtain a weak order which 1/2 in time 
and 1 in space. 

2 Preliminaries 
2.1 Functional spaces. 

It is convenient to change the notations and rewrite the unknown of p.2p as X. We thus 
consider the following stochastic partial differential equation written in the abstract form 

(2.4) dXt + AXt dt = Q^/^dWt, Xq = x € H, 0<t<T, 

where H is a Hilbert space whose inner product is denoted by (•,•) and its associated 
norm by | • |, the process {^t}tg[o,T] is an ff-valued stochastic process, A a non negative 
self-adjoint unbounded operator on H whose domain D{A) is dense in H and compactly 
embedded in H, Q a non negative symmetric operator on H and {VFt}j(=[o,T] ^ cylindrical 
Wiener process on H adapted to a given normal filtration {J't}te[o,T] i^i ^ given probability 
space (J^,J^,P). 

It is well known that there exists a sequence of nondecreasing positive real numbers {An}n>i 
together with {e„}„>i a Hilbertian basis of H such that 

Acn = Xne-n with lim A„ = +00. 

n— >+oo 
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We set for any s > 0, 

D{A^/'^) = }u = ^ UnCn e H such that ^ X^u^ < +00 \ , 



+ 00 +CXD 



and 



n=l n=l 



n>l 



It is clear that D{A^/'^) endowed with the norm \\u\\s '■= \A^/'^u\ is a Hilbert space. We 
define also D{A~^/'^) with s > as the completed space of H for the topology induced by 
the norm = Y2n>i Ki^'^'ri- this case D{A^^/'^) can be identified with the topological 

dual of D{A'I'^), i.e. the space of the linear forms on D^A^/"^) which are continuous with 
respect to the topology induced by the norm || • \\s. 

Moreover, theses spaces can be obtained by interpolation between them. Indeed, for any 
reals si < s < S2, one has the continuous embeddings C D{A^/'^) C D{A^'^/'^) and 

by Holder inequality 

(2.5) hlls < s = (1-A)si + As2, 
for any u G 

We denote by || • ||x the norm of a Banach space X. If X and Y denote two Banach 
spaces, we denote by C{X, Y) the Banach space of bounded linear operators from X into 
Y endowed with the norm = ^^Vx&x ll-^2;||y/||x||x- When X = y, we use the 

shorter notation C{X). 

If L G C{H) is a nuclear operator, Tr(L) denotes the trace of the operator L, i.e. 

TY{L) = Y^{Lei,ei) < +00. 

It is well known that the previous definition does not depend on the choice of the Hilbertian 
basis. Moreover, the following properties hold 

(2.6) Tr(LM) = Tr(ML), for any L, M € C{H), 
and 

(2.7) TV(LM) < D:{L)\\M\\^H), for any L G C+iH), M e C{H), 

where JC+{H) denotes the set non negative bounded linear operators on H. 

Hilbert-Schmidt operators play also an important role. Given two Hilbert spaces Ki, K2, 
an operator L G C{Ki,K2) is Hilbert-Schmidt if L*L is a nuclear operator on Ki or 
equivalently if LL* is nuclear on K2. We denote by C2{Ki, K2) the space of such operators. 
It is a Hilbert space for the norm 

m\c,(K„K,) = {TT{L*L)f' = {Tv{LL*) f' . 

It is classical that, given four Hilbert spaces Ki, K2, K3, K4, if L G jC2iK2, K^), M G 
7C{Ki,K2), N G C{Ks,K4) then NLM G C2{Ki,K4) and 

(2.8) \\NLM\\c,^K,,K,) < \mciKs,K,)\\L\\c,{K,,Ks)\\Mhc{K,,K,)- 



4 



See [6], appendix C, or [TO] for more details on nuclear and Hilbert-Schmidt operators. 

If X is a Banach space, we denote by Cb{H;X) the Banach space of X-valued, continuous 
and bounded functions on H. We also denote by C^{H) the space of /c-times continuously 
differentiable real valued functions on H. The first order differential of a function (p G 
C^{H) is identified with its gradient and is then considered as an element of Cb{H;H). It 
is denoted by Dip. Similarly, the second order differential of a function ip G C^{H) is seen 
as a function from H into the Banach space C{H) and is denoted by D'^ip. 



2.2 The deterministic stationary problem 

We need some classical results on the deterministic stationary version of ()2.4p . In this case, 
special attention has to be paid to the space V = Di^A^/"^) C H. It is a Hilbert space whose 
embedding into H is dense and continuous. Its inner product is denoted by ((•,•)). We 
have 

{{u,v)) = {A^/'^u,A^/'^v), for any u G y, w G y 

= (^4^,^), for any u G D[A), v £ H. 

Then by a density argument and the uniqueness of the Riesz representation (in V) we 
conclude that A is invertible from V into V' = D{A~^^'^) or from D{A) into H. We will 
set T = A-^ its inverse. It is bounded and positive on H and on V. 

For any f € H , u = Tf is by definition the unique solution of the following problem 

(2.9) uGV, {{u,v)) = {f,v), for any t; G 

Let {Vh}h>o be a family of finite dimensional subspaces of V parametrized by a small 
parameter h > 0. For any h > 0, we denote by Ph (resp. Ilh) the orthogonal projector 
from H (resp. V) onto Vh with respect to the inner product (•, •) (resp. ((•, •))). 

For any /i > 0, we denote by A^ the linear bounded operator from Vh into Vh defined by 

(2.10) {{uh,Vh)) = {AhUh,Vh) = {Auh,Vh) for any Uh G Vh, Vh G Vh- 

It is clear that Ah '■ Vh ^ Vh is also invertible. Its inverse is denoted by T/j. For any f & H, 
Uh = Thf is by definition the solution of the following problem : 

(2.11) Uh G Vh, {{uh,Vh)) = {f,Vh) = {Phf,Vh), for any Vh G Vh- 

It is also clear that Ah and Th are positive definite symmetric bounded linear operators on 
Vh- We denote by {\^h}i<i<i(h) the sequence of its nonincreasing positive eigenvalues and 
{ei,h]i<i<i{h) the associated orthonormal basis of Vh of its eigenvectors. Again, by Holder 
inequality, Ah satisfies the following interpolation inequality 

(2.12) \Aiuh\ < \A;iuh\^\A'i^Uh\^-\ Uh £Vh, s = Xsi + (1 - X)s2. 
The consequences of ()2.10p are summarized in the following Lemma. 

Lemma 2.1 Let Ah G C{Vh) defined in ^KM)- Let T and Th defined in (KW and [2J1]) . 
Then the following hold for any Wh Vh and v G V: 

(2.13) ThPh = UhT, 

(2.14) \A'/^Wh\ = lAl^^Whl 

(2.15) \Tl^^Wh\ = \T^^W\, 

(2.16) \Th^^Phv\ < \T^^^v\. 
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Proof 

Now let f £ H. We consider the two solutions u and Uh of ()2.9p and (I2.11|) . Since Vh C V, 
we can write (|2.9p with Vh S Vh- Then, substracting we get {{u — Uh,Vh)) = 0. Hence, 
Uh = II/jU the F-orthogonal projection of u onto Vh, i.e. ThPhf = ^hTf. 

Equation (12.140 follows immediately from the definition (I2.10p of Ah- We now prove ()2.16p . 
Equations (|2T3]1 and (|2JiD imply \A~^^^Phv\ = \\ThPhv\\ = \\IlhTv\\ < \\Tv\\ = \A-'^/'^v\, 
since Uh : V ^ Vh is an orthogonal projection for the inner product || • ||. 

As regards (|2.15p . on one hand ()2.16p with v = Wh & Vh C V gives the first inequalty 
\Th^Wh\ < \T'^^'^Wh\- On the other hand, by (|230]) . 

\{AhUh,Vh)\ = \{Auh,Vh)\ < lA'/'uhWA^Wl = \A'/'uh\\A^/'vh\- 

So AhUh can be considered as a continuous linear form on D{A^^'^), i.e. belongs to D(A~^/'^), 
and 

\A-y'AhUh\<\AyW\ = \Al/\h\. 
Taking Uh = A~^Wh gives \ A-^/^Wh\ < Eq. ([213]) follows. | 



Our main assumptions concerning the spaces Vh is that the corresponding linear elliptic 
problem (j2.1ip admits an 0{h^) error estimates in H and 0{h'^^^) in V for some r > 2. 
It is classical to verify that these estimates hold if we suppose that Uh satisfies for some 
constant hq > 0, 

(2.17) \UhV-v\ < Koh'\A'/'^v\, l<s<r, 

(2.18) \A^/^(nhW-w)\ < noh'''^\A''^'^w\, 1 < s' < r - 1, 

where v E D{A'/'^) and w G D{A'''/'^). 

Finite elements satisfying these conditions are for example Pk triangular elements on a 
polygonal domain or Qk rectangular finite element on a rectangular domain provided A; > 1. 
Approximation by splines can also be considered. (See [3], |3U| ) . 



2.3 The deterministic evolution problem 



We recall now some results about the spatial discretization of the solution of the determin- 
istic linear parabolic evolution equation: 

du(t) 

(2.19) -^ + Au{t) = 0, n(0)=y, 



by the finite dimensional one 

dt 



+ AhUhit) = 0, UhiO) = PhV e Vh- 



It is well known that, under our assumptions, ()2.19p defines a contraction semi-group on 
H denoted by S{t) = e"*"^ for any t > 0. Its solution can be read as u{t) = S{t)y where 
t > 0. The main properties of S{t) (contraction, regularization) are summed up below: 



(2.20) |e"*^x| < \x\, for any x e H, 
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and 

(2.21) \A'e-^^x\<C{s)t-'\x\, 

for any t > 0, s > and x € H. Such a property is based on the definition of and the 
following well known inequality 

(2.22) supx^e-*"" <C{e)t-', for any t > 0. 

x>0 



In the same manner, we denote by Sh{t) or e the semi-group on Vh such that Uh{t) = 
Sh{t)Phy, for any t > 0. 

We have various types of convergence of Uh towards u depending on the regularity of the 
initial data y. The optimal rates of convergences remain the same as in the corresponding 
stationary problems (see ()2.17p ~ (l2.18p ). The estimates are not uniform in time near t = 
since the regularization of S{t) is used to prove them. The following Lemma gives two 
classical properties needed in this article. 



Lemma 2.2 Let r > 2 be such that (2.11) and ^2.18\) hold and q, q, s, s such that 
0<s<q<r,s'>0, l<q' + s'<r — 1 and q' < 2. Then there exists constants Ki > 0, 
i = 1,2 independent on h such that for any time t > 0, one has: 

(2.23) \\Shit)Ph-smciDiA'f^),H) < ^ihH-^'i-^y^ 

(2.24) l|5h,(t)P/.-5(t)IL(,,(^.'/2),z5(Ai/2)) < 



The proof of ()2.23|) can be found in [2j (see also ^32j, Theorem 3.5, p. 45). The proof 
of (|2.24|) can be found in [21], Theorem 4.1, p. 342 (with / = 0)). In fact, we use only 
s = s' = and q' = 1 below. 



2.4 Infinite dimensional stochastic integrals 

In this section, we recall basic results on the stochastic integral with respect to the cylin- 
drical Wiener process Wt- More details can be found for instance in |6]. 

It is well known that Wt has the following expansion 

+ 00 

m = J^A(t)ei, 

i=l 

where {/3j}i>i denotes a family of real valued mutually independent Brownian motions on 
{fl,T,F,{Tt}t>o)- The sum does not converge in H and this reflects the bad regularity 
property of the cylindrical Wiener process. However, it converges a.s. and in L'P{Vt; U), 
p > 1, for any space U such that H C U with a Hilbert-Schmidt embedding, li H = L'^{0), 
O dW^ open and bonded, we can take U = H^%0), s > d/2. 

Such a Wiener process can be characterized by 

E{Wt,u){Ws,v) = mm{t,s){u,v) 
for any t,s > and u,v £ H. 
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Given any predictable operator valued function t ^{t), t £ [0,T], it is possible to define 
^{s)dW{s) in a Hilbert space if $ takes values in C2{H, K) and \\^{^)\\']:.2{h K)*^^ ^ 
CO a.s. In this case ^{s)dW{s) is a well defined random variable with values in K and 

[ ^s)dWis) = y^[ ^s)em{s). 
Jo Jo 

Moreover, if E (f^ \\^{s)\\'j;_^(^jj j^-^ds] < oo, then 



and 



We will consider below expressions of the form jQij{s)Q^/'^dW{s). These are then square 
integrable random variables in H with zero average if 



Ie(^ \ms)Q'^^\\l,iH,K)ds^ Ti {r{s)Qib(s))ds<^. 

The solution of equation ()2.4|) can be written explicitly in terms of stochastic integrals. In 
order that these are well defined, we assume throughout this paper that there exists real 
numbers a > and min(a — 1,0) < P < a such that 

oo 

(2.25) ^ = \\A-/Yc,iH) = Tr(vl-") < +oo. 

n=l 

and 

(2.26) A^Q G C{H). 

Condition p. 260 implies that Q is a bounded operator from H into D{A^). By interpolation, 
we deduce immediately that for any A G [0, 1], A-^^Q^ G C{H) and 

(2.27) \\A'(^Q>^\\c{H) < WA^Qfc^HY 



Example 2.3 // one considers the equations described in the introduction where A is the 
Laplace operator with Dirichlet boundary conditions, it is well known that Ii2.25\) holds for 
a > d/2. 

We have the following result. 

Proposition 2.4 Assume that ( fE^) . / TO^) hold and 
(2.28) l-a + p>0. 

Then there exists a unique Gaussian stochastic process which is the weak solution (in the 
PDE sense) of {2.4^ continuous in time with values in L'^(^}, H). It is given by the formula 
which holds a.s. in H: 

ft +°o / ft \ 

Xt = e-'^x + e-(*-^)V/'dW^, = e-*^x + ^^ e-(*-)^Mft(s)J Q^/^a. 
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Proof 

By Theorem 5.4 p. 121 in [6|, it is sufficient to see that the stochastic integral make sense 
in i.e. 

for any t £ [0,T]. We use ()2.8p to estimate the Hilbert-Schmidt norm: 

< c{t - s)-V2(-/3+a) 

by (^THh . (IT^ and (IT^ . The conclusion follows since -/? + a < 1. I 



3 Weak convergence of an implicit scheme. 
3.1 Setting of the problem and main result. 

In this section, we state the weak approximation result on the full discretization of (j2.4p . 

We first describe the numerical scheme. Let > 1 be an integer and {Vh}h^Q the family 
of finite element spaces introduced in Section [2l Let At = T/N and tn = n At, < n < N . 
For any h > and any integer n < N, we seek for XJ^, an approximation of Xt^, such that 
for any Vh in 14: 

(3.29) iXj^-^'-Xj:,Vh) + AtiAXi:+',VH) = {Q'^'Wt^^, - Q'/^Wt„,Vh), 
with the initial condition 

(3.30) {Xlvh) = {x,Vh), yv^eVh, 
where 

= exj^+^ + (1 - 0)x]i, 

with 

(3.31) 1/2<6'<1. 

Recall that for 9 < 1/2, the scheme is in general unstable and a CFL condition is necessary. 
Then ()3.29p - (l3.30p can be rewritten as 

(3.32) X^+i -XJi + At AhXl^' = VAtPhQ^I\''+^ , 

(3.33) Xl = Phx, 
where 

V At 

and where we recall that P/^ : if — > V/j is the ii-orthogonal projector. Hence {x"}n>o is a 
sequence of independent and identically distributed gaussian random variables. The main 
result of this paper is stated below. 
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Theorem 3.1 Let Lp S C^{H), i.e. a twice differentiable real valued functional defined 
on H whose first and second derivatives are bounded. Let a > and fi > be such that 
112. 25\) . 112. 26\) and /12. 28\) hold. Let T > 1 and {Xt}t^[o^T] be the LL-valued stochastic process 
solution of ll[2.4\ ) given by Proposition \2.4\ For any N >1, let {-^^}o<n<Af be the solution 
of the scheme k3.3i3^) - l3.'J3\) . Then there exists a constant C = C{T,Lp) > which does not 
depend on h and N such that for any 7<1 — a + /3<l, the following inequality holds 

(3.34) |E fiXj^) - E <f{XT) I < C (/i^T + At^) , 
where At = T/N < I. 

3.2 Proof of Theorem \3A\ 

The scheme (I3.32p - (l3.33p can be rewritten as 

n-l 

(3.35) XJ: = Sl^,PhX + ^tY,SlJ^-\l + eAAh)-^Phx''^\ 0<n<N, 

k=0 

where we have set for any h > and > 1: 

ShAt = {! + OAtAh)-\l - (1 - 6)AtAh). 

Step 1: We introduce discrete and semi-discrete auxihary schemes which wih be usefuU 
for the proof of Theorem 13.11 

First, for any /i > 0, let {Xh{t)}t&\Q^T] be the V/i-valued stochastic process solution of the 
following finite dimensional stochastic partial differential equation 

dXh,t + AhXh^t dt = Ph Q^'^dWt, Xh,o = Phx. 

It is straightforward to see that Xh^t can be written as 

(3.36) Xh,t = Sh{t)PhX + r Sh{t - s)PhQ^'^dWs. 

Jo 

The last stochastic integral is well defined since t ^ Ti {{Sh{t)PhQ^''^Y{Sh{t)PhQ^^'^)) is 
integrable on [0, T] . 

We introduce also the following V/j-valued stochastic process 

Yh,t = Sh{T-t)Xh,t, tG[0,T], 
which is solution of the following drift-free finite dimensional stochastic differential equation 

(3.37) dYh,t = Sh{T - t)Ph Q^'^dWt, Yhfl = Sh{T)PhX. 
Its discrete counterpart is given by 

(3.38) y," = Sj^^^^XJi, 0<n<N, 

n — 1 

= sf^^^^Phx + VAtY, sf^J-Hi + eAtAh)-^PhQ^/\^+\ 

k=0 
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Eventually, we consider a time continuous interpolation of Yj^ which is the V/j-valued {^t}t- 
adapted stochastic process Yh^t defined by 

(3.39) Yh,t = Sf^^^PhX + / V S^l^~\l + eMAhrHn{s)PhQ^'^dWs, 
where 1„ denotes the function Ij^^^^^^^. 

It is easy to see that for any t G [0, T] and n be such that t £ we have 

Step 2: Splitting of the error. 

Let now ip G C^{H). The error ¥.ip{Xl^) — E,ip{Xt) can be splitted into two terms: 

(3.40) E99(X^)-E(/.(Xt) = Eip{X^)-E^{Xh,T)+^^{Xh,T)-^^{XT) 

= A + B. 

The term A contains the error due to the time discretization and will be estimated uniformly 
with respect to h. The term B contains the spatial error. 

Step 3: Estimate of the time discretization error. 

Let us now estimate the time error uniformly with respect to h. In order to do this, we 
consider the solution u : V/i — > M of the following deterministic finite dimensional Cauchy 
problem: 

(3.41) I ^ = ((^^(^ - t)PhQ'/'rD\„{S,{T - t)PuQ"^] 

We have the following classical representation of the solution of (|3.41|) at any time t G [0, T] 
and for any y G V^: 

(3.42) Vh{T -t,y) = ¥.ip(^ + Sh{T - s)Ph Q^/^dW^^ . 
It follows easily 

(3.43) Wvhimci^H) < MlciiH), t e [0,T]. 

Now, the estimate of the time error relies mainly on the comparison of Ito formula applied 
successively to t i— > Vh{T — t,Yh^t) and t i-^ Vh{T — t,Yh^t)- First, by construction, t i— > 
Vh(T — t,Yh^t) is a martingale. Indeed, Ito formula gives 



dvh{T - t, Yh,t) = [Dvh{T - t, Yh,t), Sh{T - t)PhQ'^^dWt) . 

Therefore 

,(r - t, Yh,t) = Vh{T, Sh{T)Phx) + j'^ {Dvh{T - s, Yh,s), Sh{T - s)PhQ^/^dWs) . 



Taking t = T and the expectation implies 

(3.44) E^{Xh,T) = Vh{T,Sh{T)Phx). 
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On the contrary, t Vh{T — t, Yh^t) is not a martingale. Nevertheless, applying Ito formula 
gives, thanks to (j3.39p . 



rp TV— 1 

(3.45) E 

■^0 n=0 



dvh 
dt 



(T-t,Yh,t)dt 



ln(i)cli, 



where here and in equations (j3.46p . (j3.47p below, D'^Vh is evaluated at (T — t, Y/j t). Also 
we have set 

ThAt = {I + OMAh)-\ 
Now, plugging ()3.4ip into (I3.45|) gives: 

(3.46) + i E / ^Tr [s^^Zt-'THAtPkQ'^') DW {slir'TnAtPHQ'") 

•^0 n=0 L 



ln(t) dt. 



At last, the comparison between ()3.44p and ()3.46p leads to the following decomposition of 
the time error A 



^^{Xi^) - E99(X;,,t) = Vh{T,S';:^^^Phx)-Vh{T,Sh{T)P^,x) 
(3.47) +-Ey^ 5]Tr {sl-^r'ThAtPHQ"^y D\n {slir'ThAtPhQ' 



/2 



n=0 



Sh{T - t)PhQ^'^y D\h [Sh{T - t)PhQ^/^ 
I + 11, 



ln(t)dt. 



The term / is the pure deterministic part of the time error. Thanks to the representation 
(f3:i2D . we have 
(3.48) 



I < M\cl{H)\\Sh{T)Ph - Sh^i^tPh\\c{H)\x\- 

Thanks to (I3.3ip it is possible to bound / uniformly with respect to h. More precisely, we 
have 



(3.49) 



||(5,(iVAt)-5^%)P,||^(H) = sup 



< suple"^" -F^(z)| 

z>0 



for any T > 1 (e.g. see Theorem 1.1, p. 921 in [23j). We have used the following notation: 

, l-(l-e)z 
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Let us now see how to estimate the term II. First, using the symmetry of D'^Vh, we rewrite 
the trace term as 

+2 TV (^{Sh{T - t)PhQ^l^y D\h [sll^l-^n^tPhQ^I^ - Sh{T - t)PhQ^I^' ' 

= an{t) + bn{t). 



Let now q > and /3 > such that (p:25|) and (p:26|) hold with < 1 - a + /3 < 1. Let 
7 > and 71 > such that 0<7<7i<l-a + /3<l. 

We first estimate the term a„(t). We use dZS]), (pT]) . (f03]l and to obtain 



an{t) <\\D^Vh{T -t)\\c2 



< M 

< M 



Cl(H)\ 
Cl(H)\ 



, ^hAt ^'PhAt - Sh{T - t) ] PhQ^^'^\?c^(H) 



S-, 



N-n-1 
hAt 



ThAt-Sh{T-t)]A^ 



(l-7i)/2 
h 



P 



(3.50) 

Note that here Vh is endowed with the norm of H. Let us set 



Mn{t) 



(3.51) 



^hA^t ^PhAt -Sh{T -t)) A^ 



(l-7i)/2 



sup 

l<i<I{h) 



K,h{T-t) 



1 + eXi.hAt 



Ph 

x(l-7i)/2 



Using similar techniques as for the proof of the strong order (see e.g. [28]), we have the 
following bound, for n < — 1 



(3.52) 



Mn{t) < C7AtT/2((iV -n- l)At)-(^- 



-7l+7)/2 



where here and below C denotes a constant which depends only on 71, 7, Hy'll (72(^1^), 
and Tr(A~°). In particular these constants do not depend on h or At. The proof is post- 
poned to the appendix at the end of this article. 



We then estimate the last factor in (I3.50p . Since Vh C H and we have endowed Vh with the 
norm of H, we may write 



1^ 



(7i"l)/2 p ^1/2|| ^ II .(7l-l)/2 p ^1/2|| 

PhQ \\C2{H,Vn) ^ ll^h PhQ \\C2{H)- 



Using (12. 8|) . we deduce 



I 4(71-1)72 p ^1/2 II ^ 
\Af^ PhQ \\C2iH,Vh) 



C2{H) 



Using (121^ and (12^61) . we have 



1^ 



(7i-l)/2 



PkA-P/Yr 



C2{H) 



Him 



,(7:~l)/2p^^_;3/2g^|2 

< E.eN |A-^/V,^-/^/2e.p(i-70|p,^-/^/2e,|27i 

< E.e^N |^-^/2-'^/2e,|2(i-70|^-/5/2e.|2,, 

— Z^ieN \ 

We deduce from 1 - 71 + /? > a and (|2.27|) that 

\\A^r'^"PHA'^'Yc2m < Xl-''^'-\\A^Q\\Z^TriA-n 



13 



Plugging this and (13.520 into ()3.50p yields for n < — 1: 

(3.53) anit) < CAt^iiN -n- l)At)-^^-^'+^\ 
For n = — 1, we derive similarly, 

(3.54) ^ C {mAtAt'-'^Xv,, + \\S.{T - t)At'^^/%y^,) 

< C7(AtTi-i + (r-t)Ti-i) 

< c(r-t)Ti-i. 

Concerning we write 

bn{t) = 2Tr ((5,(T - t)4^-^^)/^4^^-^)/^P,Qi/2)* ^2^^ 

< Mc.^j,^\\S,iT - t)Af^/'\\c^v,)\\ {s^,Ar'TH,At - SniT - t)) Af^/'\\c^v,) 

Using similar argument as above, we prove 

\\S,{T - t)A^^-'''^/'\\c^y^^ < C{T - t)(^-i)/2 

and, for n < N — 1, 

II {sl'^r'THAt - Sh{T - t)) A^t''^^"\\c(y^) < CAt\{N -n- l)At)-(i-^^+^) 
so that n < N — 1: 

(3.55) hn{t) < CAt^{{N -n- l)At)-^'^-^^+^\ 
For n = — 1, we have 

hM-i{t) < ||v.|b2(^)||5,(r-t)4^-^^)/^||^(v.,) 

(3-56) x|| (r,,A, - S,{T - t))4'""^^^'^^k(V.)ll4'""^^^'^/.Q^/'llL(..,v.) 

< c{T-tyi^-^ 

We are now ready to bound // in (fTTTll . Indeed, ([TS^ . ([T^ . (fT^ imply 

(3.57) // < C Y,At^{{N -n-l)At)-^^-^'+^Hn{t)dt+ C{T-ty'~^dt 

J At Jo 

< CAt"' 

Then, plugging ([Tig]) and ([H^T]) into ([TITl) we obtain that 

(3.58) 1^1 < K4||/)(^||c,(//;H)|x|At+^r^i-^At^ 

< CAtT, 

for T > 1, At < 1. 

Step 4: Estimate of the space discretization error. 
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Let us now estimate the spatial error B. The method is essentially the same as above: we 
use the Kolmogorov equation associated to the transformed process Yj. 

We consider the following linear parabolic equation on H: 

(3.59) ^{t,x) = ^TT(^DMt,x){S{T-t)Q'/'){S{T-t)Q'/'r), t > 0, x e H, 

together with the initial condition 

v{0,x) = (p{x), X £ H, 

where f is a real-valued function of t and x £ H. We have the following representation of 
V (see e.g. [7], chapter 3) at time t £ [0, T] and any y £ H: 



(3.60) 



v(T -t,y) =Eip [y 



+ SiT - s) Q^^'^dWs^ . 



We apply the Ito formula to t i— > v{T — t,Yt) and t i-^ v{T — t,Yh^t)- We substract the 
resulting equations and obtain 



( EifiXr) - E^{Xh,T) = v{T, S{T)x) - v{T, Sh{T)Phx) 

T 



(3.61) < 



+^eI Tr (^Sh{T-t)P„Q'/^yDMT-t,Yh,t)[ShiT-t)PhQ'/^') 



dt 



--E 



Tr 



5(T-^)QV2^ DMT-t,Yh,t)(s(T-t 



dt. 



ShiT - t)Ph 

Sh{T - t)PhQ'/'' - S{T - m^f^y D^v{Sh{T - t)PhQ'/^ - S{T - t)Q'/^) 
SiT - t)Q^lA* D^v{Sh{T - t)PhQ^/^ - S{T - t)Qy^' 



The first term on the right hand side of ()3.6ip is the deterministic spatial error which can 
be bounded thanks to ([2:23]) (with s = and g = 27 < 2) and ([3:60]) . We obtain: 

(3.62) |z;(T,5(r)x)-z;(T,5;,(r)P/,x)| <Ki||(^|b^(;,)i/i27T-^|x|. 

For the second term, we use the symmetry of D^v and write 

Tr (^Sh{T-t)Pf,Qy^yD\(sUT-t)PhQ'/^) - ( S{T - t)Qy^y D\( S{T - t)Qy 

= Tr ' 
+ 2Tr 

= a + b, 

where here and below D^v is evaluated at {T — t, y^. 



Let 7 > be such that (I2.28|) holds and 71 > such that 0<7<7i<l-a + /3<l. 
Thanks to (|2.7p and (|3.60p . we get the following bounds: 



b = 2Tr (^S(T-t)DMSh{T-t)Ph-S{T-t))Q'^ 

= 2Tt(^A'^^~^-'^ A^-^^ S{T - t)D^v{Sh{T - t)Ph - S{T-t) 
< 2 \\{SUT - t)Ph - S{T - mc(H)MyAH)U'~^'S{T - t)\\c(H)\m^\\ciH)TT{A^'''^-P). 
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Then, owing to ^2^, (fOS]) (with s = and g = 27 < 2), we obtain 

(3.63) b < Ch^^'iT - t)-(i+T-Ti) 

where again C denotes a constant which depends only on 71, 7, Hy'llcz^//), ||^^Q||£(_f/) and 
Tt(A~°') but not on h or At. 

As regard a, we get first thanks to (j2.7|) and (j3.6Up : 

a < \\ip\\c2^H^TT(Q{Sh{T-t)Ph-S{T-t)y{Sh{T-t)Ph-S{T-t))') 
= McUH) \\Q'^'{Sh{T - m - S{T - t))||2^(^) 

< C\\A^'-'^y\S^{T - t)P, - SiT - t))||2(^), 

where we have used (I2.27p . If 71 + 7 > 1, we interpolate ()2.23p with = (7 + 71 — l)/7, 
s = and (l2:2il) with s' = 0, q' = 1 and get 

(3.64) \\A^'-^^y\Sh{T - t)Ph - S{T - t))||^(^) < C h?{T - t)-(i-7i+7)/2_ 

If 71 + 7 < 1, we interpolate (fT^ with g = 0, s = and (fT^ with s' = q' = 1 and get 

(3.65) \\A^^-^^'^'\Sh{T-t)Pt,-S{T-t))\\c(H)<Ch^-^^{T-t)-^^-^^\ 
We use again an interpolation argument to get 



\\A(^-^^)/\SH{T-t)PH-S{T-t))\\c(H) 

< C\\{Sh{T - t)Ph - S{T - tmi\^^\\A^/\Sh{T - t)Ph - S{T ^ ^^m'^"^^ 

< C {WA^^S^iT - mWciH) + \\Ay^S{T - t)\\a^H)f~'' 

< c(r-t)(i-7i)/2, 



thanks to ()2.2ip for A and A^ and p.l4p . A further interpolation between ()3.65p and ()3.66p 
gives 

(3.67) p(l-^^)/2(^^(^ _ _ g^rj. _ < C7/iA(l-7i)(2- _ ^)~(l-7i)(A+(l-A)/2)_ 

Taking A = 7/(1 — 71) shows that (I3.64p again holds for 71 + 7 < 1. 
We deduce 

(3.68) a < Ch^^{T - t)-^+^'^-^ . 

Plugging (j3.62p . (j3.63p and (j3.68p into (|3.6ip leads to, after time integration which is 
relevant since 1 — 71 + 7 < 1 : 

(3.69) \B\ < C h^^, 
for T > 1. 

Conclusion: Gathering (I3.58P and ()3.69p in (I3.40p ends the proof of Theorem 13.11 
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Appendix. 



We now prove the estimate (I3.52p on M„, n < N — 1, defined in (I3.5ip 
We proceed as follows: 



(3.70) Mn{t) < sup 

l<i</(h.) 



1 + eXi.At 



A 



(l-7i)/2 
i,h 
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+ sup 

l<i<I{h) 



-^i.h(T-t„+i) _ p—>^i.h{T—t) 



1 + ehhi^t 



A 



-A,,,,{T-t) / 



1 



+ sup e , X - I 

l<i</(h) V + f^^t^i,h 

= ai + 02 + 03. 

Thanks to p.49p with replaced by — n — 1 we get 



{l-7i)/2 
i,h 

(l-7i)/2 



A 



< 



K4 



A 



sup 



(l-7i)/2 
i,h 



- n - 1 i>f \ 1 + eAi,,,At 



{N-n- i)(i-7i+7)/2 I 1 + 0Ai,ftAt 



{l-7i)/2' 



< 



((A^-n-l)At)(i-7i+7)/2 



Indeed, since (1 - 71 + 7)/2 < 1, {N - n - I) > [N - n - i)(i-7i+7)/2, the same way, 
we have 

' - fin (i + eAtx,,,) J 



< 



< 



C(7) At^/2 sup ( Af^-^^+^)/2e-(^-«-i)^*^.'^ 

i>l ^ ' 

C(7,7i)At^/^ 
((iV-n - l)At)(i-^i+T)/2' 



where we have used that \t — tn+i\ < At and the inequahty |e ^ — e ^| < C-^lx — yp^^ and 
()2.22p . Eventuahy, similar computations lead to 

< C7(7,7i) At^/2(^_^)-(l-7i+7)/2 

< C(7, 71) At^((iV - n - l)At)-(i-^i+^)/2^ 
where we have used again the inequality ()2.22p . 

Gathering these three estimates in (I3.70p yields (I3.52p . for n < — 1. 
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